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Abstract. The energy-dependence of charged particle mean multiplicity and pseudorapidity density at 
midrapidity measured in nucleus-nucleus and (anti)proton-proton collisions are studied in the entire avail- 
able energy range. The study is performed using a model, which considers the multiparticle production 
process according to the dissipating energy of the participants and their types, namely a combination of 
the constituent quark picture together with Landau relativistic hydrodynamics. The model reveals inter- 
relations between the variables under study measured in nucleus- nucleus and nucleon-nucleon collisions. 
Measurements in nuclear reactions are shown to be well reproduced by the measurements in pp/pp in- 
teractions and the corresponding fits are presented. Different observations in other types of collisions are 
discussed in the framework of the proposed model. Predictions are made for measurements at the forth- 
coming LHC energies. 
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1. Soft hadron multiparticle production is one of the 
most intriguing topics in high-energy interaction studies. 
Data have been investigated in different types of interac- 
tions, ranging from lepton-lepton to nucleus-nucleus in- 
teractions, and over a large energy span, covering several 
orders of magnitude. QCD, the theory of strong interac- 
tions, has provided partonic description of many observa- 
tions. However, the problem of soft multiparticle produc- 
tion still eludes a complete understanding and remains one 
of the challenging problems in high-energy physics [1] . The 
new high-energy data from LHC provide an opportunity 
to look at the system under new conditions. Of special in- 
terest are nucleus- nucleus collisions, probing nuclear mat- 
ter at extreme conditions, where new forms of matter are 
expected to be created at very high densities and tem- 
peratures. The data available from RHIC experiments al- 
low an interesting comparison of the particle production 
mechanisms with the less complex e + e~ and pp systems. 
In this context, the global, or bulk, variables such as the 
average charged particle multiplicity and particle densi- 
ties (spectra), which are the first available experimental 
observables, are of fundamental interest [TJEI as tne y are 
sensitive to the underlying interaction dynamics. 

In this paper we consider the center-of-mass (cm.) 
energy dependence of the average multiplicity and near 
midrapidity density of charged hadrons produced in nucleus- 
nucleus and (anti)proton-proton collisions. Whereas the 



multiplicity is sensitive mostly to the fraction of energy 
being transformed into observed particles in a given reac- 
tion, the midrapidity density reflects different stages of 
the reaction. Both variables increase with the collision 
cm. energy. Recent measurements at RHIC follow the 
trends observed in e + e _ and proton-proton interactions. 
The values of both bulk variables are found [3l|U to be 



similar when comparing the measurements in e + e _ inter- 
actions at the cm. energy of y/s cc , and in most central 
("head-on") heavy- ion collisions at the nucleon-nucleon 
cm. energy 



where the measurements in the 
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latter case are normalized to the number of pairs of par- 
ticipants ("wounded" nuclcons [5]). This phenomenon is 
found to be independent of the type of colliding nucleus 
for y/iii between - 20 GeV and 200 GeV. 

Assuming a universal mechanism of hadron production 
is present in both types of interaction, and that it is driven 
only by the amount of energy involved into secondary pro- 
duction, one would expect the same value of the observ- 
ables to be obtained in proton-proton collisions when ^/s pp 
is almost equal to ^/snn- However, comparing these data 
[51[7] to the measurements from RHIC, one finds j3]0,8,9, 
ll^[TTl[T2l[T3] significantly lower values in hadron-hadron 
collisions. Furthermore, the recent RHIC data from deuteron- 
gold interactions at ^/snn = 200 GeV unambiguously 
points to the same values of the mean multiplicity as mea- 
sured in antiproton-proton collisions [HQ3]- 

To interpret these findings, we have proposed in [15] 
a phenomenological description based on the energy dis- 
sipation by colliding participants into the state formed 
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during the early stage of the collision. Particle production 
is then driven by the amount of the initial effective en- 
ergy deposited in this early phase by the relevant types 
of participants. The experimental observations referred to 
above have been shown to be well described by this model 
and further predictions have been made. In this paper, the 
new and higher-energy data are added and analyzed. 

2. In our consideration, the whole process of a colli- 
sion is treated as the expansion and the subsequent break- 
up into particles from an initial state, in which the total 
available energy is assumed to be concentrated in a small 
Lorcntz-contracted volume. There are no restrictions due 
to the conservation of quantum numbers other than en- 
ergy and momentum constraints, thus allowing a relation 
between the amount of energy deposited in the collision 
zone and the features of bulk variables in different reac- 
tions. This approach resembles the Landau phenomeno- 
logical hydrodynamical description of multiparticle pro- 
duction in relativistic particle collisions |16j . Though the 
hydrodynamical description does not match ideally the 
data on multiparticle production in the whole range of 
pscudorapidity and different particle species, it gives good 
agreement with the multiplicity measurements in such dif- 
ferent reactions as nucleus-nucleus, pp, e + e~ and i/p col- 
lisions demonstrating striking predictive power [^l fTSlfT7T 
H51ITgil201|g3P 5 ] . Recently, the Landau model prediction 
for the Gaussian pscudorapidity shape due to the longitu- 
dinal particle transport has been shown to reproduce well 
the RHIC data [T51I24] as well as [221125] the phenomenon 
known as the "limiting fragmentation" [25] . The latter was 
demonstrated to be independent of the energy and types 
of colliding objects |4l lT9ll20| . This indicates that the main 
assertions of the Landau approach are useful to estimate 
fractions of the energy dissipated into particles produced 
in different reactions, particularly in nucleus-nucleus colli- 
sions [2B] . Let us stress here that in this paper the Landau 
hydrodynamical model is considered in the frame of the 
constituent quark picture as it is described below. 

Once the collision of the two Lorentz-contracted parti- 
cles has resulted in a fully thcrmalizcd system, but before 
expansion, we assume that the production of secondary 
particles is defined by the fraction of energy of the partic- 
ipants deposited in the volume of thcrmalizcd system at 
the moment of collision. This implies that there is a dif- 
ference between results of collisions of structureless and 
composite particles: in composite particle collisions not 
all the constituents deposit their energy when they form a 
small Lorcntz-contracted volume of the thermalized initial 
state. Therefore, in nucleon-nuclcon collisions the interac- 
tions occur between single constituent, or dressed, quarks 
in accordance with the additive quark picture [27] , and the 
other quarks are considered to be spectators. Thus the en- 
ergy of the initial thermalized state which is responsible 
for the number of produced secondary particles is that 
of the interacting single quark pair. The quark spectators 
which are not part of the thermalized volume at the mo- 
ment of collision do not participate in secondary particle 
production. As a result, the leading particles [25] resulting 
from the spectator quarks carry away a significant part of 



the energy. Thus, only about 1/3 of the entire nucleon 
energy is available for particle production in pp/pp colli- 
sions. 

In heavy ion collisions, however, more than one quark 
per nucleon interacts due to the large size of the nucleus 
and to the long travel path inside the nucleus. The more 
central the nucleus-nucleus collision is, the more interac- 
tions occur and the larger is the energy available for sec- 
ondary particle production. In central nuclear collisions, 
a contribution of constituent quarks rather than partic- 
ipating nucleons seems to determine particle production 
and their distributions [23]. In the most central collisions, 
the density of matter is so high (almost saturated) that 
all three constituent quarks from each nucleon may par- 
ticipate nearly simultaneously in the collision, depositing 
their energy coherently into the thermalized collision vol- 
ume. In this case, the entire energy of the participating nu- 
cleons is available for bulk production in head-on nucleus- 
nucleus collisions. Comparing this to proton-proton colli- 
sions, where only one out of three constituent quarks from 
each proton interacts, one expects the features of the bulk 
variables per pair of participants measured in the most 
central heavy-ion interactions to be similar to those from 
proton-proton collisions but at a three times larger cm. 
energy, ^/s~^ ~ 3 ^s NN . 

Adding together the above discussed ingredients, namely 
the Landau model and the constituent quark picture, one 
finds for the ratio of the charged particle rapidity den- 
sity p(y) = {2/N- pay: t)dN c \ l /dy per participant pair at the 
midrapidity value y = in heavy- ion reaction, p(0), to the 
density /9 PP (0) in pp/pp interaction, 



2JV, 



ch 



P(0) 

Ppp(O) N^N™ 



L 



pp 
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Here, iV par t is the number of participants (N p 



v part 2 

in nucleon- nucleon interactions), A r c h and are the 

mean multiplicities in nucleus- nucleus and pp/pp interac- 



In t^— with m being the mass 
the mass of the proton, 



TP ' 



tions, respectively, and L 
of a participant, e.g. m - 
in nucleus-nucleus collisions. According to our model, we 
compute the ratio (pQ) for the rapidity density p(0) and the 
multiplicity AT c h at ^/snn and the rapidity density p pp (0) 
and the multiplicity JVSJ' at 3 ^/snn- Due to the above, we 
consider a constituent quark of mass \m v as a participant 
in pp/pp collisions, and a proton as an effectively struc- 
tureless participant in head-on nucleus-nucleus collisions. 
Then, from Eq. (fTJ) one obtains: 



p(0) = p pp (0) 



2^ch 

^partiV c P h P 



41n3 



ln(4m2/s NN ) ' 



(2) 

p(0) is thus calculated from the measured values of p pp (0) 
and the multiplicities measured in both reactions. 

Solving Eq. (JTJ for N ch at given p(0),p pp (0) and 
one finds: 



2iVch = ^pp P(0) 

part ° h Ppp(O) 



2 In 3 



N, 



In (4.5 v /s NN /m p ) 



(3) 
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Fig. 1. The pseudorapidity density of charged particles per participant pair at midrapidity as a function of cm. energy per 
nucleon, ^/snn, measured in central nucleus-nucleus (A A) collisions and calculated from pp/pp data using Eq. ([2}. The AA 
data are from: the AuAu measurements at RHIC by BRAHMS [22], PHENIX PMTJ, PHOBOS [5P^P5PfJl[57] . and STAR 
[3811391140] experiments; the values recalculated in [TD] from the measurements at CERN SPS by CERES/NA45 [32] and NA49 
[55] experiments, at Fermilab AGS by E802 and E917 experiments [3T], and at GSI by FOPI Collab. [3D]; the PHOBOS data 
on CuCu collisions from [36U37II42] , The solid symbols show the values obtained from Eq. (J2| using the following data on 
midrapidity densities measured in non-single diffractive collisions: pp data from ALICE 23] and CMS 03] experiments at LHC 
and from pp collisions by UA5 Collab. at CERN SPS [51135] and ISR = 53 GeV), by CDF Collab. at Fermilab 0; and in 

the following inelastic collisions: pp data by UA5 Collaboration and pp data from the LHC ALICE [43] experiment and from 
the ISR [JB] and bubble chamber [491150] experiments, the latter as recalculated in [pj. The solid line shows the linear- log fit, 
—0.33 + 0.38 ln(sNN), to the AA data with the parameters and errors obtained using a combination of the data from the RHIC 
and SPS experiments. The shaded area shows l-cr error band to the fitted parameters. The circled stars show the heavy-ion 
predictions for LHC A A collisions (open stars) and from the expected LHC pp collisions (solid stars), both calculated from the 
fit. 



From these one calculates the expected values for nucleus- 
nucleus central collisions, since they are related to pp/pp 
interactions through the constituent quarks picture of in- 
teractions combined with the Landau energy dissipation 
hydrodynamics. Note that the different sets of observables 
used in Eqs. (|2J and ([3]), do not guarantee that both equa- 
tions describe the data well, even if one of them does. 

3. In this section, we consider the available data on 
the mean multiplicity and midrapidity density measured 
in nucleus-nucleus collisions. 

Figured] shows the nuclear data on the pseudorapidity 
density, p(0), per participant pair, measured in head-on 



nucleus- nucleus collisions (open symbols), function 
of v^nn in the energy range ^/snn = 2 — 200 GeV. The 
measurements are from experiments at GSI [30] , AGS [3T] , 
CERN SPS [3133], and RHIC [SUinillllEZIEllESlISMSZI 
38,35]. The data clearly satisfies a linear- log relationship. 
Fitting, using a weighted combination of the data from 
the RHIC and SPS experiments, we find: 

p(0) = (-0.327 ± 0.026) + (0.381 ± 0.021) ln(s NN ). 

Similar observations have been made in [lTjiriTl[351r4T] . Such 
a behaviour seems [20] to be well justified by the participant- 
driven picture of the formation of the rapidity distribution 
and the iV part scaling. In addition, as the nuclear data 
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Fig. 2. The charged particle mean multiplicity per participant pair as a function of the nucleon-nucleon cm. energy, ^/snn, 
measured in the most central nucleus- nucleus (AA) collisions and calculated using pp/pp data from Eq. (J3|. The solid symbols 
show the multiplicity values from the AA data as measured by PHOBOS Collab. at RHIC in AuAu 3,4,23,34,35 and CuCu 
[36] collisions, by NA49 Collab. at CERN SPS [53] and by E895 Collab. at AGS [54] (see also [3]). The low-energy HADES 
measurements at GSI are taken from [55] . The open symbols show the values obtained from Eq. for the AA mean multiplicity 
based on: pp collisions at FNAL by E735 Collab. [5rJll57] . at CERN by UA5 Collab. at = 546 GeV [45] and v /s^ = 200 

and 900 GeV [5J]; pp collisions at CERN-ISR 05], and the ALICE, CMS and CDF multiplicities, calculated using the E735 fit, 
3.102 Sp P 178 [57], and the data from bubble chamber experiments [501159] the latter having been compiled and analysed in |60] . 
The solid line shows the 2nd-order log-polynomial fit —0.35 + 0.241n(sNN) + 0.241n 2 (sNN) to the AA data. The shaded area 
shows l-<7 error band to the fitted parameters. The dashed line shows the power-law fit: —7.32 + 5.92 s^f 1 / 4 . The stars give the 
heavy-ion predictions for the LHC AA collisions (solid symbols) and those expected from LHC pp collisions (open symbols) 
based on the 2nd-order log-polynomial fit. 



show, the process of bulk production is largely character- 
ized by iV part , or soft scaling, rather than by the num- 
ber of binary collisions, N co \\, the latter being attributed 
to hard scaling |1 01123 H M]l^5][n][5^]. Recently, this feature 
has been confirmed by observing the iV part scaling in the 
ratio of the charged-particle yields measured at ^/snn = 
200 GeV and 62.4 GeV, as (p T ) increases from 0.25 GeV 
to ~ 4 GeV [37]. The increase of the ratio at large (px) is 
believed to appear because of the harder spectra at higher 
y / SNN, so as to follow JVcoii scaling and not to scale with 
Np alt . Interestingly, the -/Vp art scaling is observed to be the 
same for AuAu and CuCu interactions. 

Figure [5] shows the ^/sNN-depcndcnce of the nucleus- 
nucleus data of the mean multiplicity iV c h per participant 
pair from the most central nucleus- nucleus collisions (solid 
symbols) in the same ^/snn range as above. The data are 



taken from the measurements by the PHOBOS experi- 
ment at RHIC [aSMlSlElESJES], NA49 at CERN SPS [S3] 
and by E895 at AGS [51] . The mean multiplicities are seen 
to increase as a 2nd-order logarithmic polynomial, and our 
fit gives: 

AT ch /(0.5JV part ) = (-0.35 ±0.28) 

+ (0.24 ± 0.16) ln(s NN ) + (0.24 ± 0.03) ln 2 (s NN ), 

which well reproduces the energy-dependence of the mea- 
surements. This trend holds even for relatively low-energies, 
as demonstrated by the recent data by HADES [55; dis- 
played in Fig. [5] Such an ln 2 (sNN) behaviour is expected 
[201I2T] to appear naturally as a combination of the fol- 
lowing features: the Landau model Gaussian shape of the 
pscudorapidity distribution, the logarithmic increase of 
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the midrapidity density with the cm. energy, and the lim- 
iting fragmentation. 

4. The midrapidity density and the multiplicity en- 
ergy dependencies are analyzed in the framework of the 
model developed here, i.e. applying Eqs. @ and © to 
calculate these variables in nucleus-nucleus interactions 
based on nuclcon-nuclcon data. 

Usingthe pp/pp data on p pp (0) from CERN [6"II35I 
l^lMM^Hpl and Fermilab [TUHEI], the data on multiplic- 
ities iV c p h p in pp/pp interactions [1^111111551155] > and N <& 
from Fig. (3] along with the fact that the transformation 
factor from rapidity to pseudorapidity does not influence 
the above ratio, we calculate, according to Eq. the 
in = a/Sdd/3. These calculations arc shown 



p(0) at 



by solid symbols in Fig. [T] One can see that the calcu- 
lated p(0) values are in a very good agreement with the 
measured p(0)'s as well as with the obtained log- fit in the 
whole measured ^/snn range. 



a preferable fit to the data [351I57] , is almost as good 
as the log 2 (sNN) polynomial approximation. However, it 
is evident that, after inclusion of the new higher energy 
data, the log 2 (sNN) function is more preferable. 



5. From the above, we conclude that using all the 
world available measurements on the mean multiplicity 
and the midrapidity density in the data of the nucleon- 
nuclcon and central nucleus-nucleus collisions, a clear in- 
terrealtion between the two types of the data is obtained. 
This can be attributed to the universality of the multi- 
particle production process over almost three orders of 
magnitude of y' snn- Under this assumption, predictions 
for the LHC energies can be made. 

Using the fits, shown in Figs. [T] and [5J and solving 
Eq. (J!} for the midrapidity density p PP (0) with from 
the high-energy fit [ST], the expected p PP (0) values for pp 
collisions at LHC are found to be about 5.8, 6.4, and 6.9 
at ^/sT™ = 7, 10 and 14 TeV, respectively, within 5% to 



The agreement is more intriping as one adds the hadronic 10% uncertainties . From the fit to the mi drapidity den 



data at cm. energies higher than the top RHIC energy 
namely the densities for (anti)proton-proton interactions 
at v /s P 7 = 1.8 TeV from Tevatron [7] and ^/s^ = 2.36 TeV 
43 .43] from the LHC. In Fig. [T] we compare the cal- 
culations using Eq. ([2]) on these TeV-energy data with 
our linear- log fit. One sees the nucleus- nucleus p(0)'s at 
^/snn = 600 and ~ 800 GeV, calculated from the high- 
est energy hadron data, agree well with the AA-fit. This 
observation justifies the above conclusion from the lower- 
energy midrapidity density dependence on the types of 
participants. Note that Eq. ([2]) shows the relevance of 
the Landau hydrodynamical energy-dissipation L-factor 
which was required in order to correctly estimate the midra- 
pidity density. 

Addressing now Eq. ([3]), we calculate the participant- 
pair-normalized mean multiplicity N^/ (0.5A" part ) for nucleus 



sities p(0), shown in Fig. [T] the p(0) values, expected for 
PbPb collisions at LHC energies ^/snn, corresponding to 



the above 



are found to be about 5.7, 6.0, and 6.2 



at y/s^ = 2.76, 3.94 and 5.52 TeV, respectively. The p(0) 
predictions are shown in Fig. [TJ by circled solid stars for 
PbPb interactions and by circled open stars for those from 

7/3 when cal- 



thc pp expectations at LHC at 
culated according to our model. 

Comparing our predictions for p pp (0) to the predic- 
tions of other models and Monte Carlo tunes [STllol] . we 
find that our values are in the range of the midrapidity 
density values predicted there. The values we find here 
are also consistent with those from the CDF p pp (0) fit [7J 
and from a similar, but higher-energy, CMS fit [33] . The 
values of p(0) we find for LHC heavy-ion collisions are also 
well in the range of the expectations by different models 



i m JST ^ r ° m ^ C PP / p Pp^rf SUr ! m f m m$m- The p(0) values obtained at ^ a 5.5 TeV is 
IBKESKflEM of P PP (0) and ACT and the corre- ^ ^ that obtained fay PHENIX V from the fit to the 

spending p(0) data from F^.ffl The calculated ^ch/(0.5^ pa rt nuclear data and b PH0B0S from their e xtrapo- 

lation of their AuAu data to PbPb collisions at LHC [20] . 

Using the log 2 snn fit to the mean multiplicity shown 
in Fig. [3] one finds the average multiplicity, N c h/ (0.5Af part ) , 
in PbPb collisions to be about 64, 70 and 73 at y / snn = 
2.76,3.94 and 5.52 TeV, respectively, with 10% to 15% 
uncertainties. The in pp collisions at LHC are ex- 
pected to be about 73, 82 and 93 at y/Sp~p~ = 7, 10 and 
14 TeV within about 10% uncertainties, and are the same 
as one finds from the multiplicity high-energy power-law 
fit [57]. The N ch /(0.5N paxt ) predictions for PbPb collisions 
are shown in Fig. [2] by solid stars, while those expected 
from the LHC pp collisions at ^/snn = v/Spp/3 in the 
framework of our model, are shown by open stars there. 

The N c h/(0.5N pp ) we obtained are consistent with ex- 
trapolations from the present experimental measurements 
|20j . Our prediction at v / snn — 5.5 TeV is comparable 
with the estimate from the pseudorapidity density spec- 
tra by PHENIX fKfllTT] . Similar to the p pp (0) values, the 
A^ are within the range of the predictions by different 
models and Monte-Carlo tunes [571I6T] . 



values are shown in Fig. [5] (open symbols) along with the 
available nuclear data. One can see that the calculated 
N c h/ (fi.5N part ) values reproduce well the log 2 (sNN) fit ob- 
tained here and follow the nuclear data points for ^/snn = 
2 GeV to about 200 GeV. 

The calculations for the new higher energy data namely 
those for -^snn ^ 200 GeV, are again of a special interest. 
These calculations help to check the proposed description 
over a larger energy range, the possible fit and result in 
firmer predictions. Indeed, from Fig. [5] one can see that 
for ^/snn ^ 200 GeV, it is quite difficult to distinguish be- 
tween the two fits: the power-law fit oc sJ N , known to be 



For the LHC data we exclude some results to avoid any de- 
pendence on the experimental cuts and to have consistent data 
sets. We do not use the 900 GeV data by ATLAS [47] which 
applies a high pr-threshold, and the inelastic (high-energy) 
data by ALICE [48] where specific selection cut is used. This 
exclusion does not influence the consideration here. 

2 For the ALICE, CMS and CDF multiplicities the E735 
power-law fit = 3.102 s° p 178 [57] is used. 
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6. Let us now dwell on some corollaries of the model 
proposed here and discuss the results in view of other ob- 
servations. 

From our consideration it follows that, at the same 
a/snn, the mean multiplicities as well as the midrapidity 
densities, normalized to the number of participants, would 
give the similar values when measured in central symmet- 
ric nucleus-nucleus collisions of different colliding nuclei, 
they arc largely driven by the initial energy deposited by 
the participants at early stage of collisions. Indeed, as seen 
from Figs. Q] and [2j this effect has been already observed at 
SPS energies and now is confirmed by the RHIC measure- 
ments at ,/snn of about 50 GeV to 200 GeV. The same 
values for both the observables are obtained [20,41,42,64 
in Au-Au and Cu-Cu data, as shown in Figs. [1] and [2J 
Note that this effect has also been observed for the whole 
pscudorapidity region [B5]. Notice also the similarity in 
the above-discussed iV part dependence in low-pr vs. high 
Pt ranges [37] . 

An interesting issue to be addressed in the framework 
of the model is to consider asymmetric collisions, such as 
nucleon-nucleus (pA/dA) ones. In such type of interac- 
tions, the bulk variables studied here, being measured at 
given y/fiNN are expected to have the same values as those 
in pp/pp collisions at ^/s pp ~ ^/snn- Indeed, assuming 
an incident proton in p-nucleus collisions interacts in the 
same way it would interact in pp collision, the secondary 
particles in the reaction are assumed to be created out of 
the energy deposited by the interaction of a single pair of 
constituent quarks, one from the proton and another one 
from a nuclcon in the interacting nucleus. This, in its turn, 
implies that the mean multiplicity and the midrapidity 
density are expected to be independent of the centrality 
of nucleus-induced collisions, V par t (within uncertainties 
due to intranuclear effects, e.g. Fermi motion). These ex- 
pectation are shown to be well confirmed in the RHIC 
data on dcuteron-gold interactions at ^/snn = 200 GeV. 
Moreover, the effect obtained at RHIC is shown [4,14, 
l3"6] to be true also for hadron-nucleus collisions at lower 
y/iii" w 10-20 GeV. 

In this study the Landau hydrodynamic model is used, 
leading to a good description of the data from different 
reactions |^[l^[T7lP^[19l[2n|[22][23l[24] . However, this is a 
1+1 model and therefore does not take into account the 
transverse expansion of the system which can be stud- 
ied, for example, via the transverse energy, another im- 
portant bulk observable. Considering the measurements 
from SPS to RHIC of the transverse energy midrapidity 
density, one finds that this variable scales with the num- 
ber of constituent quarks a way similar to that of charged 
particle and photon midrapidity densities |66j . Further- 
more, as measured at RHIC, the ratio of the transverse 
energy midrapidity density to that of multiplicity density 
is observed [TU] to be independent of the number of par- 
ticipants, and the transverse energy loss is found [57] to 
be independent of the type of colliding nuclei if the same 
number of participating nucleons is taken. The observa- 
tions indicate scalings of a similar nature to those con- 
sidered here for multiplicities and midrapidity densities. 



This seems also to reflect the fact that the inclusion of 
the transverse expansion in the Landau model does not 
change the scaling of the observables under study [TTllrjg] . 
Currently, the model attracts high interest and has under- 
gone a generalization, see e.g. [55]. 

As we have shown, the constituent quark is a key com- 
ponent of a correct description of "soft" particle observ- 
ables in particle and nuclear collisions from a few GeV up 
to highest LHC energies, so that constituent quarks have 
to be taken as the interacting particles, see e.g. [70]. This 
already has support from studies of the multiplicity dis- 
tributions in heavy-ion-induccd interactions [29,64,65,71, 
172"] . The constituent quark picture has been exploited to 
reasonably model the heavy-ion pseudorapidity and trans- 
verse energy data [73] . The elliptic flow characteristics are 
also observed to scale when the constituent quark frame 
is taken into account [741175] . 

The similarities of pp/pp and A A interactions observed 
here for the two basic variables suggest that the system is 
formed at early stage as a superposition of contributions 
from the constituent quarks. The particle multiplicities 
seem then to be derived by the total energy of partic- 
ipants, available in the Lorentz-contracted volume. The 
importance of the very early stage of collision for soft par- 
ticle production has been already discussed elsewhere [65l 
I7TII75"1[76] . Due to the proportionality of the multiplicity 
to the entropy [16] . the multiplicity scaling, observed for 
different variables measured at RHIC, has been suggested 
to be connected to the total produced entropy [77] . 

Considering this and recalling the above-mentioned 
similarity in heavy-ion and e + e _ collisions observed, one 
would expect the same model to be valid in matching the 
mean multiplicity and the midrapidity values in heavy-ion 
and e + e~ data. Indeed, as we have shown in [T5], both 
variables follow the same energy dependence within the 
framework of our picture, as soon as one considers that 
the structureless electron and positron deposit their total 
energy into the Lorentz-contracted volume similar to nu- 
cleons in central nuclear collisions. From this, the factor 
1/3 applied to the pp/pp energy scale is expected to re- 
sult in a good match between the e + e _ and pp/pp data 
on multiplicity and midrapidity densities as shown in |15| . 
This solves the problem with the energy-scaling factor of 
1/2 used in OS], where the yfs pp /2 shift is shown to pro- 
vide a reasonable description of the average multiplicity 
cm. energy dependence but not of the midrapidity den- 
sity when comparing heavy- ion/c + c _ data to those from 
pp/pp collisions. We recall that the energy-scaling factor 
1/3, has already been shown in [751175"] to give good agree- 
ment of the pp mean multiplicity data relative to those 
from c + e _ annihilation, for a review see pQ. It is remark- 
able that the 3NLO perturbative QCD [80] fit to e+e" 
data [SI] describes the pp/pp multiplicity data providing 
the inelasticity is set to ~ 0.35 [ST], favouring the effec- 
tive 1/3 cm. energy in multihadron production in pp/pp 
reactions. 

As already mentioned, the average multiplicity is de- 
fined mostly by the fraction of the cm. energy transformed 
into observed particles, so that, after the energy shift in 
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y/s P p is applied, the pp/pp data reproduce reasonably well 
the e + e _ multiplicity data. For the mid-rapidity density, 
the subsequent system development has to be taken into 
account and is shown well described by the Landau hy- 
drodynamics picture providing the cm. energy is scaled 
according to the contribution of the participants [i"5ll79] . 

7. In summary, we analyse the average multiplic- 
ity and midrapidity density data in pp/pp and in cen- 
tral nuclear interactions as a function of the cm. energy 
per nucleon over the whole available range of the interac- 
tion cm. energies, including the highest energy LHC data 
y/Spp = 2.36 TeV in pp/pp collisions, and the highest en- 
ergy RHIC data at ^/snn = 200 GcV. Within the frame- 
work of constituent quarks, we develop a model which 
interrelates these two variables measured in the two types 
of interactions, assuming one quark of each nucleon par- 
ticipates in pp/pp collision while all three quarks (i.e. a 
complete nucleon) participate in a head-on heavy ion col- 
lision. We consider these participants to form the initial 
zone of a collision which then develops in hydrodynamic 
framework, the Landau relativistic hydrodynamic model 
in our case. In this approach soft hadron production is 
determined at the very early stage of the collision. Af- 
ter appropriately taking into account the contributions of 
the participants, which requires an energy-scaling factor 
of 1/3 in pp/pp measurements, the average multiplicity 
and mid-rapidity density in nucleon-nucleon and nucleus- 
nucleus interactions are found to have a similar cm. en- 
ergy dependence. The midrapidity density is found to obey 
a linear-log fit on ^/snn, while the multiplicity data fol- 
lows a second-order log-polynomial increase with ^/snn- 
A clear preference of the data on the multiplicity to follow 
the log 2 (s) behaviour is observed compared to the power- 
law, the two dependencies being indistinguishable up to 
^/snn of about 300 GeV, or y/s pp ~ 1 TeV. Assuming no 
changes in the multihadron production processes with in- 
creasing energy i/Sp~p~ of the LHC to 7 TeV, 10 TeV and 
14 TeV, and looking forward to the heavy-ion data at the 
corresponding of 2.76 TeV, 3.94 TeV and 5.52 TeV, 

we estimate the multiplicities and midrapidity densities 
for the forthcoming data, using the obtained energy de- 
pendencies. 

We are grateful to David Plane for his help during preparation 
of the manuscript. 
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